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, The exterior solution for an arbitrary charged, massive source, is studied as a 

static deviation from the Reissner-Nordstr0m metric. This is reduced to two 
, coupled ordinary differential equations for the gravitational and electrostatic 

m . 

, potential functions. The homogeneous equations are explicitly solved in the 

ON , 

particular case q 2 = m 2 , obtaining a multipole expansion with radial hyper- 

ON 



cr 



geometric dependence for both potentials. In the limiting case of a neutral 
source, the equations are shown to coincide with recent results by Bondi and 
Rindler. 



X 

H 

03 1 I. INTRODUCTION 



A series of well-known and elegant results in general relativity (see references in 
show that the multipole moments of an isolated, stationary source determine in a complete 
and unique manner the gravitational field on its surrounding 3-space. This has been also 
extended to the case of charged sources, both for static || and stationary |§ cases. Recently 
however, Bondi and Rindler stated the convenience of obtaining precise expressions for the 
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multipole expansion of a fully- relativistic isolated source by means of a static perturbation 
defined over the Schwarzschild metric. In this way, comparison between newtonian gravity 
and general relativity can be achieved in quantitative terms. In the same spirit, this work 
presents an extension a la Bondi and Rindler to the Einstein-Maxwell case. Thus, we 
consider a bounded and static, but otherwise arbitrary distribution of charged matter. The 
assumed static character ensures that the only conserved quantities for such a system are 
the total gravitational mass-energy m and the net charge q, as evaluated by an observer 
situated well away from the source. It is sensible to expect that such an observer would 
find a space-time not markedly different from the Reissner-Nordstr0m one. As Bondi and 
Rindler would say, "the Reissner-Nordstr0m metric gives bones to the space" . We therefore 
propose a description of the exterior geometry as a small static deviation from the Reissner- 
Nordstr0m solution. 

This paper is organized as follows: in section [T1J the relevant Einstein-Maxwell equations 
are presented, and the constraints imposed by the field equations on the perturbation terms 
are found. In section [Lll| , the physical implications of the results are interpreted. The 
particular cases q 2 = m 2 and q = are discussed in detail. 



In curvature coordinates (t,r,9,<p), and taking advantage of the static character of the 
system so as to eliminate the mixed time-space components of the metric tensor, and to 
diagonalize its pure spatial part by means of a coordinate transformation, we are allowed to 
adopt the following general ansatz for the line element: 



II. THE EINSTEIN-MAXWELL FIELD EQUATIONS 





where 




(2) 
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is the usual Reissner-Nordstr0m metric function, and a, (3, 7 and 5 are small functions of 
the spatial coordinates. On the other hand, the vector potential for the electric field will be 
chosen as 

*a = (- ~ *) 5° a , (3) 



with $ defined again as a small function of the spatial coordinates. It is obvious from ([!]) 
and (0) that the five-parameter set of perturbations a, /3, 7, 5 and $ represent all possible 
non-spherical terms for the physical fields, and that their vanishing brings us back to the 
spherically symmetric Reissner- Nordstr0m solution. This set of functions, defined in a 
somewhat similar manner to [J], is also completely equivalent to the time-independent part 
of Chandrasekhar's polar perturbations M , when considered as first-order quantities. Their 
behaviour is governed by the electrovacuum Einstein-Maxwell field equations, 

Rab = 8nT a b 

= 2F ac F\ + ^g ab F cd F cd , (4) 
V b F ab = 0, 

where F a b = 2V[ a Aq is the electromagnetic Maxwell- Faraday tensor. Since we are interested 
in the far-field solutions, we can neglect squares and products of the perturbations, as well 
as of their derivatives, so that the non-vanishing Einstein-Maxwell equations become a set 
of eight partial differential equations, which we transcribe below: 
(R tt = SnTu): 
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Resolution 



Since all static and axisymmetric solutions to the Einstein-Maxwell equations belong 
to Weyl's electrovac class 0, it is convenient to define the five perturbation functions in 
such a way as to avoid the repetition of these known results. On the other hand, it is also 
convenient to exclude the possibility that those perturbations arise simply as a consequence 
of an inappropriate choice of coordinate origin, in which case even the spherically symmet- 
ric solution would possess such non-zero metric functions. Both exclusions may be easily 
achieved in the selected coordinate system, by means of an expansion in spherical harmonics 
in the usual way: 



oo n 



a(r,9,p)=Y: E °ha(r)Yl(6,<p), (13) 

n=2 l=—n 
(1^0) 

and other similar expansions for the remaining functions /3, 7, 5 and $. It is important to 
notice that the chosen set of indexes n, I explicitly imposes the above mentioned restrictions, 
for n = represents the spherically symmetric terms, which are already fully present in the 
background Reissner-Nordstr0m curvature; the I = set corresponds to the axisymmetric 
multipoles; and finally, a non- vanishing n — 1, I — ±1 term simply means a shifted choice 
of coordinate origin. 

Keeping this in mind, we notice that the field equations (JTTJ) , ( PH| ) and (|9|) impose the 
following relations between the perturbation functions 



P = —a, 
dr 

5 = 7; (14) 



2o ^ da p' 
pr z or p 



whose inclusion in the remaining field equations show that (H) and (|6|), as well as ((?]) and @ 
become identical, leaving us with just three independent equations. Using the expansions 
(Pf), and adopting the notation a (r), 7 (r) and $ (r) to represent the set of radial coefficients 
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Oinh Ini and all of them with the same index n, these equations reduce to a system of 
two ordinary differential equations on a and namely: 
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plus an algebraic relation for 7 
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In this way, a simultaneous solution 



of ([TB]), would determine, via fljif) , (p!4|) 



a (r) ; $ (r 

and the expansions (|T3|) , the Einstein-Maxwell fields (|l]) and (^). 

The so obtained field equations show the convenience of the assumed linearization in 
a Reissner-Nordstr0m background, since the absence of cross-terms has eliminated the l- 
dependence in ( |15"D for each (2n)-set of multipoles with the same index n, thus recovering 
the well-known fact, that it is always possible to restrict ourselves to an axisimmetric de- 
pendence, completely determined by the index n. 

In the limiting case of null electric field, q = $ = 0, our set of independent equations ([15]) 
reduces only to one differential equation on the gravitational perturbation a (r), namely 



d 2 a 2 (r — m) da 4m 2 + n (n + 1) (r — 2m) r „ 
-T7T H ; z — c— — -o ol = 0. 



(17) 



dr 2 r (r — 2m) dr r 2 (r — 2m)' 

The simple transformation from our radial curvature coordinate r to the radial isotropic 
coordinate p through a new dimensionless radial parameter x, 



x 



m 



m / \2 
— frad 1 + y/xj \/x, 



(18) 



converts ( p.5a| ) into 

,d 2 a 1 



, 1 X ) X 



da 



, a H — (1 — 3x) (1 — x) x , 
dx 2 2 y M ; 



4x + -n (n + 1) (1 — x) z 



a = 0, 



(19) 



for {0 < x < 1} corresponding to the region {2m < r < 00}, which is identical to the one 
obtained in Obviously, the form of (0), which is an homogeneous differential equation 



with three regular singularities, ensures that it is possible to obtain exact solutions directly 
from it, but their demonstrated equivalence with the Bondi-Rindler solutions doesn't make 
it worth while. 



III. ANALYSIS OF THE EQUATIONS 

We would like to emphasize that the electric and gravitational fields are completely 
determined by (|15|) , which are a coupled system of linear differential equations. The general 
relativistic coupling between these fields is obvious, since both equations possess on their 
right hand side an inhomogeneous term. Therefore, the general solution of (|l^) can be 
split into an homogeneous part for each function a and $, plus the inhomogeneous, coupled 
terms. 

Furthermore, ( |15aj ) and ( |15b| ) are elliptic equations, and hence they define a well-posed 



boundary problem, whose boundary conditions must be imposed as physically meaningful 
constraints. In addition, their elliptic character offers an extra advantage, for the required 
boundary conditions will have to be satisfied not only by the entire solutions a (r) ; $ (r) , 
but also by their homogeneous parts a (r) ; $ (r) , and this can be simply checked out by 
means of the simpler set of homogeneous equations, given by the left hand sides of ([151). 

The most obvious boundary condition was already assumed in the definition of the 
perturbations, for these must vanish at large distances from the distribution in order to 
recover the Reissner-Nordstr0m solution. In this way, it is straightforward to verify that 
both a and $ are determined by two arbitrary integration constants, and therefore, for both 
fields, each set of harmonics possesses the same multiplicity of solutions as the classical 
study of independent gravitational and electrostatic potentials. Likewise, in each case only 
one of the two possible solutions will vanish as they approach infinity. On the other hand, 
the second boundary condition should be imposed by the junction with a convenient interior 
solution; but its arbitrary character would make it very difficult to define in a tractable way. 
Fortunately, the metric (|I|) inherits a physical feature from the chosen background curvature, 
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at least for the general cases q 2 < m 2 , since the analytic extension of the Reissner-Nordstr0m 
solution describes a black hole, with an event horizon situated at r + = m + ^rri 1 — q 2 , and 
its "no-hair" property brings a suitable boundary condition. This can be easily shown 
analytically, and without invalidating the generality of the above remark, if we restrict 



ourselves to the simpler case q = m . The homogeneous parts of (131) can be written, in 



m 



terms of a convenient dimensionless radial parameter follows: 

r 



x 2 (1 — x) 2 
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a = 0, 
$o = 0. 



(20a) 
(20b) 



In this way, the r — > oo boundary condition is now located at the finite value x = 0, and 
the event horizon for the maximal charged black hole at x — 1, so that the outside space is 
covered by the region {0 < x < 1}. As both equations now possess three regular singularities 
at x = 0, 1, oo, the former are granted to have solutions with an hypergeo metric dependence 
given by some functions h (x), defined as follows: 



a = x Pl (1 - x) qi h x (x) , $ = x P2 (1 - xf 2 h 2 (x) 



(21) 



whose insertion in ( P0| ) brings up the following pairs of solutions for the exponents p, q, 
suitably expressed in terms of an index j = n + \: 



n+1 



Pl,P2 



Qi 



n 



--±^•2 + 2, q 2 = +- ± - 4. 



(22) 



It is then obvious that only the solutions with p equal to n + 1 will satisfy the asymptotic 
vanishing condition. On the other hand, the g's with positive root must be discarded, for 
they would lead to vanishing perturbations at the event horizon x = 1, showing a spherically 
symmetric configuration at the finite radial value r = m; this situation would be in flagrant 
contradiction with Birkhoff's theorem. In this way, we are led to a unique set of choices; 
with these, ( |2~0D can be cast into the canonical form of standard hypergeometric equations 



x (1 — x 
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(23a) 
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whose integrals are given by the hypergeometric functions 



h\ (x) — H (ax, bx, 2n + 2; x) , h 2 (x) = H (a 2 , b 2 , 2n + 2; x) 



(24) 



a 2 + b 2 = 2(n + l-v^^4) 



(25) 



with the a's and 6's defined according to 

ai + 6! = 2 (n + 1 - Vj^Ti) 
a x b x = 2 (n + 1) (j - V/FT2) | a 2 6 2 = 2 (n + 1) (j - VF^) 

Although in both cases a and b are complex numbers, the regularity of h\ and h 2 is completely 
ensured by their hypergeometric definition 

(a) n (b) n x r 



rr/ . , , a&x a(a + 1)6(6 + 1) ~ 
g(a,i,,c;x) = l + - Ij+ ^ V + '-- = S 



n! 



(26) 



n=0 ( C )' 

since all a and 6 solutions of (|25|) are complex conjugates, the functions in (J2fj) are always 
real- valued. It is then obvious that h% and h 2 possess a regular behaviour in the [0; 1] 
range; moreover, their hypergeometric dependence ensures also that they are monotonously 
increasing functions in that same range, as can be easily seen by keeping in mind their 
vanishing at x — 0, and evaluating them and their first derivatives at x — 1. So, the radial 
dependence of a and $ will be mostly dominated by the x p (1 — x) q terms in For 
x — > 0, fl22l ) shows that the asymptotic behaviour will be determined by the x n+1 term, thus 
inducing a radial dependence given by 

1 



(27) 



which corresponds to the classical coulombian decreasing rate; this also shows that the cou 
pling terms between both fields, i.e. those who are non-coulombian contributions, asymp 
totically vanish with greater powers of the radial coordinate, and only the harmonic 
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part of the fields will be perceptible. On the other side, both a and $ diverge in the x — ► 1 
limit, since their radial dependence will be dominated by the (1 — x) q terms, and the g's 
defined in (|22| ) are all negative for the permitted values of n. This divergent boundary 
condition at x = 1 is in agreement with the "no-hair" theorem for charged black holes, since 
it ensures that there can be no bounded static perturbations of the analytically extended 
electrovac solutions on the event horizon. It also shows that in the x ~ 1 strong-field region, 
the coupling between both fields dominates over their coulombian parts, given by (f27|). The 
homogeneous functions «o; $o so obtained also offer an algorithmic way to compute the 
complete inhomogeneous field equations, since these functions can now be substituted into 



. It- 



the right hand sides of ([jjj) in order to find some first inhomogeneous functions 
erating this procedure leads to more accurate inhomogeneous functions, which, when added 
to the homogeneous ones, will bring exact solutions to the field equations. Anyway, it must 
be born in mind that these functions might not be exact solutions of the field equations out 
of the weak-field condition, where their smallness condition has to be relaxed. Nevertheless, 
as discussed in 0, the ellipticity of (|1"5| ) ensures that no solutions other than these can exist 
outside of a sphere which contains all the sources. 
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